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A presentation about concepts rather than techniques 


Richardson’ s Forecast 


Table 5: Six-hour Changes in Pressure Thickness 

Richardftou's Values: No Initialization 


Layer 


Horizontal 

Convergence 

Vertical 

Convergence 

I 

48. 5 

65.9 

-17.4 

II 

28.4 

-23.7 

52.1 

in 

25. 3 

47.6 

-22.3 

IV 

22. 3 

7.5 

14.8 

V 

20.8 

48.0 

-27.2 

Sum 

145.4 

145.4 

0.0 


Lynch, Peter, 1994a: Richardson's Marvellous Forecast. Proceedings of the 
International Symposium on the Life Cycles of Extratropical Cyclones, Bergen, 
Norway, 27 June--1 July, 1994, 38--48. 


Geostrophic Adjustment 


0 hour 



Daley 1992 


1 hour 



3 hours 



2 hours 



6 hours 



CCMIGuR INTERVAL: «. u -» REPRESEMS OS U/S 


Figure 6.3 (a) Geostrophic adjustment of initial geopotential perturbation, (b-e) Solutions at 
1. 2, 3. and 6 hours. Contoured field is geopotential, and wind arrows indicate 
speed and direction of windfield. (After Barwell and Bromley, 1988) 


Filtered Equations 


lst-order balance 

e.g., quasi-geostrophic equations 

2 nd -order balance 

e.g., nonlinear balance equation and 
quasi-geostrophic omega equation 


Dynamical Initialization 





Figube 1 . — Schematic representation of iteration methods for 
initialization with the primitive forecast equations. 


Nitta and Hovermale 1969 MWR 


Normal Modes of the Linearized PE 


Linearize P. E. about u = 0, v = 0, T = T r , p s = = 0: 


Solutions: 
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Daley 1991; Errico 1986 NCAR/TN-277 
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Linear Normal-Mode Initialization 

g(t=0) = 0 


Structures 
of two 

normal modes 


Temperton and Williamson 1979 




Nonlinear 

Normal 

Mode 

Initialization 


dg/dt (t=0) = 0 


(b) 



Figure 10.1 Time evolution of surface pressure during a 24 hour model integration for (a) 
linear and (b) nonlinear normal mode initialization. Solid curves, uninitialized; 
dashed curves, initialized. (After Williamson and Temperton, Mon. Wea. Rev. 109 : 
745, 1981. The American Meteorological Society.) 


Machenhauer 1997 Cont. Atmos. Phys.; Baer and Tribbia 1997 MWR 


Why is the extra-tropical atmosphere quasi-balanced? 

Chamey: 1955 Tellus (a paraphrase) 

The observed extra-tropical motions are dominantly quasi-balanced 
because: 

1 . The principal atmospheric forcing is large scale and long period. 

2. The quasi-balanced motion must be relatively stable with respect 
to gravity- wave perturbations, (by inference; also see Errico 1981) 

3. Dissipation must be sufficient to remove what energy is otherwise 
leaked into gravity waves, (added by R. Errico ) 

Lorenz 1980 JAS Atmospheric dynamics lies on a slow-manifold. 


Gravity Waves as 

Forced and Damped harmonic Oscillators 

The amplitude of a gravity wave structure is governed by 
an equation of the form: 

dg_ 
dt 


iXg + F(t) - vg 


Consider harmonic forcing F(t) — F(0) exp(— i/it). Then 
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QG Theory and NNMI 


In the extra- tropics, the NNMI balance condition dg/dt =0 is equivalent to 

(1) The nonlinear balance equation relating mass and vorticity fields, 
with some additional small terms; 

(2) The QG-omega equation defining the wind divergence, 
with some additional small terms; 

(3) Solved with the constraint that a form of linearized potential vorticity 
is specified; 

(4) And applied only to large vertical but small horizontal scales for 
which the resonant frequency is large. 

The choice of constraint and scale selectivity matter! ! 


Harmonic Dial for External m=4 Mode, Period=3.7h 



With NNMI 



Errico 1997 J Japan Met Soc 


Harmonic Dials from a Climate Simulation 


External Mode P=3.7h 


Internal Mode P=11.6h 



Errico 1997 J Japan Met Soc 


3 


= —iujjCj + A(r,r) + B(r,g ) + C{g,g ) + D 



NATURAL FREQUENCY |cr| (s' 1 ) 

Errico 1984, 1990 MWR; Errico et al. 1988 MWR 


Why does balance matter in data assimilation? 


1 . Unrealistic initial imbalances will create unrealistic forecasts 

2. Unrealistic imbalances can be accentuated through moist diabatic processes. 

3. Large initial imbalances will tend to create less accurate backgrounds 

4. Balance can be exploited to relate u, v, T, ps (esp. in extra-tropics) 

5. Errors in balanced initial conditions will tend to create balanced background 
errors, so the error statistics should reflect that; i.e., background errors of 

u, v, T, ps tend to be correlated, esp. in extra- tropics. 


Consistency between analysis and initialization 




Fig. 2. The 500-mb height field on 14 March 1990 (a) as analyzed by NOGAPS, (b) analysis increments, (c) initialization increments, 
and (d) the sum of analysis and initialization increments. Contour interval is 60 m in (a) and 10 m in (b)-(d). Zero contours are omitted; 
negative contours are dashed; and labels in (a) are dekameters. 


Errico, Rosmond, Goerss 1993 MWR 


Global mean squared divergence tendency 
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Vertical Modes: 0 to 3 



Structures of 3 largest 
scale vertical normal 
modes 
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Lynch and Huang 1992 MWR 
Fillion et al 1995 Tellus 


IAU (Incremental Analysis Update) 


For n = — N , . . . , N — 1 
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Bloom et al. 1996 MWR 


Localization in Ensemble DA 


Example response to a single observation of § 
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Same, with 
Schur product 
Localization 





Kepert 2009 QJRMS 


Lessons Learned 


1. There are many ways to balance models, each with varying degrees 
of success. 

2. Most balance schemes have some undesirable consequences. 

3. Balance should not be applied everywhere, at all scale, in the same 
way, to the same degree. 

4. Balance should be considered when performing an analysis. 


5. Details matter. 


Common Misconceptions About Balance 


F: Small scales are not balanced. 

T: Balance depends on both vertical and horizontal scales. 
T : Deep modes are likely balanced even on the mesoscale. 

F: Atmospheric fields are on a “slow manifold.” 

T: Some atmospheric forcing has short time scales. 

T: In realistic models, freely propagating gravity waves 
are present to some degree. 


Common Misconceptions about Initialization 


F: Initialization is inappropriate when gravity waves are 
important. 

T: It is necessary when gravity waves may affect fore- 
casts. 

T: It removes waves which are not really there. 

T: It is unnecessary when gravity waves are unimpor- 
tant. 


Techniques may come and go, 
but fundamentals remain (almost) forever. 

(Unless, of course, they are neglected.) 
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